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Abstract 

Two main algorithmic approaches are known for making Hironaka's 
proof of resolution of singularities in characteristic zero constructive. Their 
main differences are the use of different notions of transforms during the 
resolution process and the different use of exceptional divisors in the de- 
scent in ambient dimension. In this article, we focus on the first difference. 
Only the approach using the weak transform has up to now been success- 
fully used in implementations, because the other one requires an explicit 
stratification by the Ifilbert-Samuel function at each step of the algorithm 
which is highly impractical due to the high complexity of the computa- 
tion of such a stratification. In this article, a (hybrid-type) algorithmic 
approach is proposed which allows the use of the strict transform without 
the full impact of the complexity of the stratification by the Hilbert- 
Samuel function at each step of the desingularization process. This new 
approach is not intended to always be superior to the previously imple- 
mented one, instead it has its strengths precisely at the weak point of the 
other one and is thus a candidate to be joined with it by an appropriate 
heuristic. 

1 Introduction 

Existence and construction of a desingularisation has been one of the cen- 
tral questions in algebraic geometry since the end of the 19th century. 
In characteristic zero, it was proved in Hironaka's groundbraking work 
[7] in 1964, in which he also introduced standard bases w.r.t. a local 
ordering among other tools, whereas the case of positive characteristic 
is still open. Nevertheless, the interest in resolution of singularities in 
characteristic zero did not end at that time. Instead, the main interest 
only shifted toward the quest for a better, more constructive understand- 
ing of Hironaka's non-constructive proof in which the key is the choice 
of appropriate centers for the blowing ups which provide the desingular- 
ization. This development led to two main algorithmic approaches to the 
proof: one due to Bierstone and Milman (see e.g. [2]) using a stratification 
by the Hilbert-Samuel-function at the beginning of each choice-of-center 
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step and the strict transform as the corresponding notion of transform. 
The other one, due to Villamayor with many contributions and simpli- 
fications by others over the years (see e.g. [T], [1], [S]), does not need 
the Hilbert-Samuel function in the choice-of-center step, but pays for it 
by using the weak transform which (compared to the strict transform) 
picks up extra components lying inside the exceptional divisor at each 
blowing upQ This latter approach has been used in the two currently 
existing implementations, mainly because its invariant is more accessible 
to practical calculations (see e.g. [3] and [^). These implementations, 
on the other hand, have brought certain questions and even conjectures, 
which had been treated purely theoretically up to that point, into the 
reach of computer experiments: among these e.g. explicit computation 
of the topological zeta-function in the quest for a counter example to the 
monodromy conjecture, treatment of singularities in singular learning ma- 
chines and in hidden Markov models in algebraic statistics or systematic 
study of multiplier ideals for cases beyond plane curves (see e.g. [11] , 
[1U|). These applications also showed that there are classes of examples, 
where the current implementations are rather far from choosing the opti- 
mal one among different possible sequences of centers. Additionally, the 
orders/degrees of the higher order generators of the ideal in question tend 
to grow much faster for the weak transform than for the strict transform. 
As the order, Villamayor's main invariant, is blind to the order of higher 
order generators, this has positive influence on the algorithmic proof, but 
very negative one on the efficiency of the implemented algorithm. 
Therefore it is a legitimate question to ask why the other approach has 
not been implemented up to now. Here the stratification by the Hilbert- 
Samuel function or more precisely the task of finding the stratum of maxi- 
mal Hilbert-Samuel function turned out to be the crucial issue. Computer 
experiments (see e.g. [S]) showed that even in very careful implementa- 
tions this step is far from being sufficiently efficient to be used in every 
choice-of-center step, because it involves a parametric standard basis cal- 
culation with as many parameters as there are variables in the basering. 
Additionally, the internal differences in the descent in dimension, which 
were already mentioned in a footnote above, turn out to be without signifi- 
cant impact on the overall efficiency, because the fewer blowing ups needed 
by Bierstone and Milman have to be paid for by the algorithmically more 
involved use of a special coordinate system and logarithmic derivatives. 
Thus a direct implementation of the Bierstone- Milman approach does not 
promise any improvement in speed or size as compared to the implemented 
algorithms. If, however, the calculation of the Hilbert-Samuel stratum can 
be avoided in a significant number of steps, this approach could again be 
an option. In this article, we follow this line of thought and provide an 
outline of an algorithm which does not compute the Hilbert-Samuel stra- 

^This difference is not the only one between the two approaches. They also differ in their 
descent in dimension of the ambient space. For Villamayor's approach we can basically use any 
local coordinate system at a given point and use the usual notion of derivatives. Bierstone and 
Milman, on the other hand, pay special attention to a choice of the local coordinates, taking 
into account the exceptional divisors passing through the given point, and use logarithmic 
derivatives to ensure that they can factor out larger powers of the exceptional divisors after 
the descent in dimension. 
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turn in each step, but instead computes auxilliary ideals which provide 
the information, whether the Hilbert-Samuel function of the original ideal 
dropped, in terms of their order. The geometric idea behind this is to 
modify the original ideal by adding extra components which only empha- 
size the contribution by certain generators of the ideal without having 
negative impact on the descending induction on dimension of the ambient 
space, which is the key to Hironakas proof and all constructive variants 
of it. 

In the section [S] we give the definition of our modified coefficient ideal 
and then apply it to the problem of resolution of singularities in the fol- 
lowing section. It is important to observe at this point, that the modified 
coefficient ideal is just another way of stating that we pass through an 
intermediate auxilliary ideal before entering the first descent in dimen- 
sion. Section 5 and 6 are then of a more practical nature, the first one 
discussing the computational issues of this new approach and the last one 
illustrating the construction on two simple explicit examples. 
As the ideas leading to this article developed in several steps over time 
and underwent more than one metamorphosis before taking the shape 
in which they are presented here, I am indepted to many people for 
fruitful discussions, among them Bernard Teissier, Vincent Cossart, Her- 
wig Hauser, Edward Bierstone, Roushdi Bahloul, Gerhard Pfister, Rocio 
Blanco and Patrick Popescu-Pampu. I'd also like to thank the Institut de 
Mathematiques de Jussieu for their hospitality, at which I was a visitor 
when my first ideas in this direction evolved. 

2 Basic Definitions and Notations 

In order to fix notation, we would like to first state the basic defintions and 
some selected properties of the invariant controlling the choice of centers. 
A section of just 1 or 2 pages can obviously not suffice to even give a brief 
introduction into the intricacies of resolution of singularities, instead we 
would like to point to more thorough discussions in section 4.2 of [5] from 
the practical point of view and in [4^ embedded in a detailed treatement 
of the resolution process. 

In the general setting for these definitions, is a smooth equidimensional 
scheme over an algebraically closed field K of characteristic zero and X C 
W a subscheme thereof. For the purpose of quickly fixing notation in this 
section, we immediately focus on one affine chart U with coordinate ring 
R and denote the maximal ideal at j: € ?7 by mx\ xi, . . . ,Xd are a local 
system of parameters of i? at a;. 

The order of an ideal I — {gi, . . . , g,-) C -R at a point x £ U is defined as 

ordx{I) := max{m G N| / C m™}; 
the locus of order > 2 of / can be computed as the vanishing locus of 

A(J) = {{g.\l <^<s}U{p^\l<^<r,l<j< d}), 
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the locus of order > c as the one of A'^ ^(^)- 

The Hilbert-Samuel function of R/I at x is defined as 

HS^ : N — > N 

s I — > lengthiR/m"^^), 

comparison is performed lexicographically. Locally at a point it can be 
computed by determining a standard basis at this point w.r.t. a local 
degree ordering and subsequent combinatorics on its leading ideal. De- 
termining its locus of maximal value, however, is an expensive parametric 
Groebner basis calculation involving as many parameters as there are vari- 
ables. 

Having fixed notation for these invariants, we can now state the general 
structure of the invariant which controls the choice of center in the reso- 
lution process: 

[ord or HS, n; ord, n; ord, n; . . .) 

where n denotes a count of certain exceptional divisors which is not going 
to have any impact on the considerations in this article. The semicolon 
in the invariant denotes the key step in the invariant, Hironaka's descent 
in dimension of the ambient space by means of a hypersurface of maximal 
contact. Hypersurfaces of maximal contact are defined by order 1 elements 
p£ ^max^ord^{i)-i^j'^^ which Satisfy certain normal crossing conditions; 
choosing a hypersurface of maximal contact can be interpreted as locally 
choosing a main variable. The order directly after the semicolon then 
denotes the order of an auxilliary ideal, the coefficient ideaff], arising in the 
descent in dimension by suitably collecting the coefficients of the powers 
of the main variable. More precisely, let Z = V{z) be a hypersurface of 
maximal contact for I at x, then the coefficient ideal of / w.r.t. Z can be 
computed as 

Coeffz(I) = 

where It is the ideal generated by all polynomials which appear as coef- 
ficients of z'' in some element of I. 

Having stated the basic notions necessary to fix notation for discussing 
algorithmic choice of center, we also need to briefly consider notions of 
transforms under blowing ups: Let tt : U — > [/ be a blowing up at a 
center C with exceptional divisor E. Then the total transform of X C U 
(defined by the ideal Ix) under the blowing up tt is given by IxOjj. The 
strict and weak transform can then be computed as 

Ix, strict ~ [IxOjj : I'e) and Ix,weak = {IxOfj : Ie) 

where b is the largest integer such that 1% ■ (IxO^ — IxOfj. 

■^More precisely, the order of the non-monomial part of the coefficient ideal is taken after 
decomposing into a product of ideals consisting of a monomial part, i.e. a product of powers 
of the exceptional divisors, and a non-monomial part. 
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3 Defining a Modified Coefficient Ideal 



In this section, we define our modified coefiicient ideal at a point w W 
in a very special situation first and subsequently study whether we can 
always pass from the general case to this particular setting. In the special 
situation, order reduction for the appropriately marked modified coeffi- 
cient ideal is equivalent to a drop in the Hilbert-Samuel function of the 
original ideal under strict transform. 

To this end, let j/i, . . . , y„ be a local system of parameters for Ow,m and 
let Tx,w = (/i, ■ ■ ■ , /fe) be subject to the following conditions: 

(1) /i, . . . , /fc is a reduced standard basis oiTx,w with respect to a local 
degree ordering such that j/i > 3/2 > ■ ■ ■ > y-a- We assume these to 
be numbered such that di < d2 < ■ ■ ■ < where di = ordw(fi). 

(2) There are integers 1 < ei < ■ ■ ■ < < n such that for each 1 < i < fc 
all V{yi), 1 < I < Ci, are hypersurfaces of maximal contact for the 
ideal 

J, = {f'"'h, . . . , y"""' /,_!, /„..., M a*") G Ir,m <r<i-l) 

where the numbers are maximal with this property and Ir,i is the 
set of all multi-indices for which X]j=i '^j'^' — 'rnax{Q,ds — dr} for 
all 1 < s < i and la'"^' | = di — d,-. 

To be able to reference coefficients of each fs £ Ji w.r.t. monomials in 
the variables i/i, . . . , j/e; separately, we write fs — '^^p '^^ Vii)^ ■ 

Definition 1 The modified coefficient ideal ofIx,w ts then defined as the 
usual coeffictent ideal of the ideal with respect to Z = V{yi, . . . , j/efc). 
More explicitly, 

Coeff-nXx,^) = £ if^ 
j=o 

where Ij = (a^'"'"' | Yllti Pi < j - dk + di). 

As will be discussed in detail in section |4l this modified coefficient 
ideal is in no way intended to be used in all descents in dimension of the 
ambient space in the computation of the value of the invariant. Instead, it 
only replaces the usual coefficient ideal in the descent of highest ambient 
dimension, accommodating for the use of the strict transform. In all fur- 
ther descents, we use again the usual coefficient ideal, which for the time 
being is the one in Villamayor's approach, although the approach could 
easily be changed to use the one of Bierstone and Milman. 

Remark 2 The use of a standard basis in the above definition is not 
necessary for the coefficient ideal itself: as all elements of a standard basis 
can be expressed as linear combinations of the original generators over 
the base ring, all contributions to the coefficient ideal already originate 
from these. The fact that we do not require identiy, but only < in the 
condition concerning the Ij again emphasizes this point of view. In the 
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subsequent section, however, the standard basis property will be used when 
considering the effects of blowing up on the coefficient ideal, because this 
allows a significantly simpler treatment of the strict transform ofTx,w 

The crucial issue about the above definition is its possible dependence 
on a particular choice of the local system of parameters. This, however, 
does not prevent its use in a resolution invariant as the following lemma 
observes: 

Lemma 3 Different choices of the yi, . . . ,yn (subject to conditions (1) 
and (2)) affect neither the order of the modified coefficient ideal nor the 
orders of the subsequent (usual) coefficient ideals. 

This, obviously, implies that the order of the non-monomial part is 
also unaffected by those different choices, as the order of the monomial 
part is not sensitive to it. 

Proof: First we observe that for Jk itself different choices of the hyper- 
surfaces of maximal contact do not affect the orders of the subsequent 
coefficient ideals, since we are using precisely the usual Villamayor-style 
resolution invariant. What remains to be proved is that the values of 
the invariant coincide for any two ideals Jk and arising from the orig- 
inal ideal Tx,w as described above w.r.t. two different local systems of 
parameter^ xi, . . . ,Xn and yi, . . . ,yn. To this end, we first consider hy- 
persurfaces of maximal contact for Ji which give rise to hypersurfaces of 
maximal contact for Ji+i and then proceed by comparing the (usual) co- 
efficient ideal of Jk and J'k with respect to suitably chosen flags. 

Step 1: descent in ambient dimension for Jk 

Let xi, . . . , j:„ be the regular system of parameters chosen in the construc- 
tion of Jk and let V{y) be a hypersurface of maximal contact for Ji for 
an arbitrary 1 < i < k (satisfying y G A'*'~^(Ji) and y ^ mod nt^„). 
Recalling that A'''~^(Ji) is generated by the generators of Ji and all their 
partial derivatives up to the (di — l)st ones, we now consider the analo- 
gously constructed set of generators of A'''+i~^(Ji+i). For those elements 
originating from standard basis elements fa with s > d, we find again all 
generators we already had plus additionally higher derivatives thereof; for 
the elements originating from an /s, s < e^, we now check the correspond- 
ing property by explicit calculation: 

Let g = be an arbitrary element of A'''~^(Ji) and let 7 = (71 -f 

di+i - di, 72, . . . ,7n)- Then 



where the summand for / = d^+i — di is precisely the desired y up to a 
non-zero constant factor. But y\ itself is an element of A''*~^(Ji) which 
implies that it can be written as a finite linear combination 



•^These two systems of parameters are of course both subject to the conditions (1) and (2). 
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where the at are constants, \rit\ < di — 1. As the appearing derivatives 
of the fi^ for it > e, are in A'''+i~^(Jf+i) by construction, we may as- 
sume without loss of generality that it < e; for all t. Replacing the /jj by 
(^a-^i-d )i. yi'^^ d, j,^^ ^j^^ ^ corresponding linear combi- 

nation yields: 



= (I-I-C2) -j/i, 

where C2 is a positive constant. Hence yi G A'**+i^"'^(Ji+i) which in turn 
implies by the above considerations that indeed g G A'''+i^^(Jf+i) as was 
to be proved. 

This shows that indeed A''^~^{Ji) is contained in A'*'+i^"'(Ji+i) which 
implies that j/s G A'''+'^^^ ( Ji+i), 1 < s < and hence proves the claim 
that V{ys) is a hypersurface of maximal contact for Ji+i. 

Step 2: Hi- Xi € mw,w 

Coming back to our original problem of comparing the coefficient ideals 
of Jfc and Jfc, we split our considerations into two parts. In this first one, 
we assume that gi := yi — xi G mfi/ j^,, and use that x\, . . . , Xe^. give rise 
to hypersurfaces of ma^ximal contact for and x\ — pi , ... , Xe^ — ge^ do 
so for Jk ■ A coordinate change replacing Xi by Xi + ffi for all 1 < i < 
transforms Jfc into and Xi—gi into Xi, thus passing from one case to the 
other. Therefore the appearing orders of all subsequent coefficient ideals 
coincide. 

Step 3: general case 

By applying the construction of step 2 as a preparation step for one of the 
two systems of parameters and as a postprocessing step for the other in 

the general case, we may restrict our considerations to a linear change of 
coordinates such that both systems of parameters satisfy the conditions 
(1) and (2). In this case, we do not even need to worry about the orders, as 
the coefficients of the monomials of a given degree are only transformed 
into linear combinations thereof still in the same degree (according to 
appropriate matrices with constant entries having full rank). Hence the 
corresponding orders are unchanged. 

q.e.d. 

After studying the effects of different choices of the system of parame- 
ters on the resulting resolution invariant, we now discuss how to pass from 
the general situation to the special situation by a variant of the standard 
basis algorithm: 

Let Tx,w be as above and let di := ordw{^x,w) be its order. Then we 

know that Jx,m contains at least one generator of order di at w. We fur- 
ther know that ei := dimKi^''^^^ {1x,w)/va^Y^^) > 0. Thus we can choose 
j/i, . . . , G mw,m giving rise to a basis of this finite dimensional vector 
space and extend this to some local system of parameters. Expressing the 
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generators of w.r.t. the new system of parameters, we enter the stan- 
dard basis calculation choosing a local degree lexicographical ordering on 
the set of monomials in the chosen system of parameters. In the standard 
basis algorithm, we only treat s-polynomials of pairs with original leading 
monomials in degree di , postponing all calculations in higher degrees. Af- 
ter appropriate renumbering and interreduction of the generators, we may 
assume that fi, . . . , fi are precisely the standard basis elements of order 
di = ■■■ = d, and satisfy LM[fx) > LM{f2) > ■■■ > LM{fi), where each 
of these leading monomials is some product of powers of j/i , . . . , i/ei . 
We now set ej := ei for all 1 < j < i- From now on, we proceed degree 
by degree through the standard basis algorithm, adding new generators 
fj to our evolving standard basis as necessary, defining the corresponding 
orders dj accordingly and setting ej := ei, until a leading monomial ap- 
pears which is not a product of the yi, - ■ ■ ,yei- We do not yet add this 
polynomial fs+i to the evolving standard basis, because we first need to 
take care of adjusting our local system of parameters in the current degree 
di. To this end, we form a second ideal J by dropping the elements of 
our evolving standard basis from the set of generators of the ideal and 
subsequently adding all products y"fj, where fj is one of the previously 
dropped elements of the evolving standard basis and |a| = di — dj. By 
construction, we now have fis+i ciimK(A'*i + '°'~^(J)/mfj/^) > = ei. 
As before, we choose j/ei+i, ■ . • , J/e^ such that yi, . . . ,ys^ give rise to a 
basis of this vector space, and extend it to a local system of parameters. 
From here on, we proceed as before degree by degree through the standard 
basis calculation extending the standard basis, defining new di and ei as 
needed and passing to a new local system of parameters according to the 
above construction when appropriate. 

Remark 4 In an affine setting, a reduced standard basis can be computed 
in family along each given connected component of a given stratum w.r.t. 
the Hubert- Samuel function. This fact will be very useful for the compu- 
tational point of view, which we will consider later on in section\^ 

4 The Modified Coefficient Ideal and the 
Resolution Process 

From the construction introduced in the previous section, it is obvious 
that a center determined in this way will always be contained in the locus 
of maximal value of the Hilbert- Samuel function of Tx ■ To study the use 
of the modified coefficient ideal in the resolution process, we would like 
to compare the effects of a blowing up in such a center to our original Jx 
and to the auxilliary ideals Jk- 

Lemma 5 The maximal value of the Hilbert- Samuel function of the strict 
transform of Ix,w is smaller than the original one at w if and only if 
the maximal order of the weak transform of Jk is strictly less than dk = 
ordyj{Jk). 
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Proof: If the maximal value of the Hilbcrt-Samuel function of the strict 
transform has decreased under the current blowing up, let us consider an 
arbitrary point wi in the preimage of w under the blowing up. There is 
at least one element, say g, of the reduced standard basis whose order 
has dropped, because the use of the reduced standard basis allows us to 
compute the strict transform by considering the strict transforms of the 
elements of the standard basis. But this implies that the order of the 
strict transform of each generator y"g of Jk is strictly less than dk and 
hence the order of the weak transform of Jk can no longer be dk- 
To prove the converse, let us now assume that the maximal order of the 
weak transform of Jk is no longer and let us fix an arbitrary point 
wi in the preimage of w under the blowing up. There is at least one 
generator of Jk, say whose strict transform has order less than dk at 
wi. A priori two situations may have occured: ordwi(hstrict) < ord^ih), 
which directly implies a drop in the maximal value of the Hilbert-Samuel 
function under this blowing up, or ord-wi {{y") strict) < \oi\. In this second 
case, we can obviously find one yj whose strict transform has order zero. 
But this yj was chosen as a local equation of a hypersurface of maximal 
contact for Tx.w and hence occurs as a factor of at least one term of lowest 
order in at least one standard basis element, say g of Tx,w But this 
implies ordmi [g strict) < ordw{g) which again corresponds to a decrease of 
the Hilbert-Samuel function as claimed. 

q.e.d. 

Lemma 6 // the maximal value of the Hilbert-Samuel function is un- 
changed under blowing up, i.e. if HSw{Ix) = HSwi{{Ix)strict), the weak 
transform of Jk coincides with the newly constructed Jk of the strict trans- 
form 0flx,w 

Proof: As the construction of the Jk uses a reduced standard basis as 

its starting point und the strict transform may be computed on the level 
of the strict transforms of the elements of a reduced standard basis, we 
immediately see that the di and the d in the construction are unchanged 
and the hypersurfaces of maximal contact may be chosen to be the strict 
transforms of the previous ones. By direct calculation of the weak trans- 
form of Jk, we can now check that the two ideals indeed coincide. 

q.e.d. 

Using the previous lemmata, we can now state a modified resolution 
invariant, which implicitly uses the Hilbert-Samuel function as the first 
building block, but only requires the explicit computation of its maximal 
value upon each drop of maximal order of an auxilliary ideal: 

(ord(Jfe),ni; ord(Coeff"''"'),n2;ord{Coeff^),n3; ...) 

in contrast to a Villamayor-style invariant 

{ord(Ix,n,),ni;ord{Coeff^),n2; ...) 

or a Bierstone-Milman-style invariant 

{HS{Ix,w), ni; ord{Coeff'^'^),n2; ...), 
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where in each case the rii denote counting of certain exceptional divisors 
and Coeff should be seen as a symbolic notation for a coefficient ideal of 
Villamayor and Bierstone-Milman respectively. 

From a theoretical point of view, the only advantage of this new resolution 
invariant lies in the fact that we may now treat the use of the strict 
transform in the framework of order reduction. From the practical point 
of view, on the other hand, this modified invariant allows us to exploit 
order reduction of higher order generators of the ideal speeding up the 
resolution process in a large class of examples without trading the fewer 
blowing ups for the huge computations necessary to determine the Ifilbert- 
Samuel stratum in each step. 

5 Computational Aspects 

For explicit calculations, it is most convenient to pass to an affine open 
covering of our smooth equidimensional scheme W and consider Ix{U) C 
Ow(U) on each of the affine open sets. For simplicity of presentation, 
we will assume from now on that Ow{U) = ^[j/i , . • . , J/n]l3 To define 
the modified coefficient ideal, which allows the use of the strict instead 
of the weak transform, our construction will proceed by iteration of two 
steps: We first need to determine (the next variety in) a suitable flag in 
W and (the next parameter in) a regular system of paramters xi, . . . ,Xn 
for W subordinate to this flag. With respect to a local degree order- 
ing on the set of monomials Mon{xi, . . . , Xn), we determine (an inter- 
mediate result up to the current degree of) a reduced standard basis of 
fi,---,fk £ 1x{U) := I C K[xi, . . . , x„]. This standard basis provides 
the complete information about the Hilbert- Samuel function through com- 
binatorial reasoning on the staircase; the coefficients of the fi w.r.t. to 
an appropriate subset of {xi, . . . ,Xi} C {a;i, . . . ,Xn} will subsequently be 
used to compute the modified coefficient ideal. 

More precisely, the construction of the flag and the standard basis can 
be stated as follows: 

StepO: Initialization 

Let di be the order of the ideal /. Create a copy lorig of / for later use. 

Stepl: Find appropriate hypersurfaces 

As the order of 7 is di, the order of A'^^~^{I) is one and any order- 
1-element thereof (subject to the appropriate normal crossing conditions) 
provides a hypersurface of maximal contact for /. If necessary cover the 

''The same construction can also be carried out in the case of Ow{U) = K[yi, . . . , yn]/ J 
for some ideal J. In that case, however, it may be necessary to pass to yet another open cover 
of U such that the local system of parameters for W at each point of the fixed smaller open 
set can be induced by the same set of elements. Even given such a set inducing the local 
systems of parameters on the whole open set, the subsequent computations still become by 
far more technical, obstructing the view to the heart of the considerations. 
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afRne chart by finitely many open sets Vi such that on each of these we 
may use the same hypersurface at all points of the locus of maximal or- 
der V{A''^~^{I)). This hypersurface gives rise to our first coordinate yi. 
Similarl£], we can proceed to determine i/2, • • • , J/ei , if there are further 
generators of order 1 of A''i~^(7). 

As we will use these new yi as the first elements of the set inducing at each 
point of the Hilbert-Samuel stratum a local regular system of parameters, 
it is convenient to use them as new variables replacing existent ones, be it 
directly, by passing to an open covering, by finite extension of our ground 
field or an appropriate combination of these methods. 

Step2 : do SB in 'degree' di 

Now we express (at least) the order di terms of the generators of I in 
terms of these new yi. We then proceed through the standard basis algo- 
rithm in this degree di in the sense that we form spolys of all pairs arising 
from generators in this degree and reduce them w.r.t. these generators 
until they are either reduced or themselves of order at least di -|- 1. 

Step3: proceed to next degree 

We note for all generators of lorig which appeared as degree di genera- 
tors of / that the respective degree is di. Then we replace all order di 
generators of / by all possible products of one such generator with one 
of the already chosen yj . Hence the new ideal I arising in this way has 
order di -I- 1, with which we can return to steps 1 and 2 with the following 
modifications: the previously chosen yj stay unchanged and the next d^ 
is defined, if additional yj arise - in addition to newly formed spolys also 
older ones are reduced, whereever necessary. Step 3 can then be applied 
using the new dj and we return to step 1 again unless we have just marked 
the last generator of lorig- 

If all generators of lorig are marked, then we have found all contribu- 
tions to the modified coefficient ideal, because all further standard basis 
element only provide coefficients which are combinations of coefficients 
already provided by lower order generators of lorig- This is sufficient for 
our purposes and we can hence stop at this point. (Note that this is not 
the complete standard basis computation and hence we cannot detect the 
whole Hilbert-Samuel function from it, but only the first entries up to 
degree dj,.) 

Step4: form modified coefficient ideal (cf. section |3]) 

On each chart which arose in the construction, we have now determined 
Ck new smooth hypersurfaces V{yi), . . . , V{ye^) and add n — ek further 
ones such that it gives rise to a regular system of parameters at each point 
of V{{yi, . - . jUe^))- With respect to these yj we can then determine the 
modified coefficient ideal. 

Note that for this coefficient ideal, we can then proceed as in the original 
algorithm of Villamayor. 

^To this end, reduce A'*i~^(/) w.r.t. the new j/i and subsequently consider its order again. 
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6 Examples 



The first example is mostly intended to illustrate what the respective 
ideals look like and how they are transformed. 

To see that the modified approach can really contribute to an improvement 

of the performance of the resolution algorithm for certain classes of ideals, 
we subsequently state a very simple explicit example for which we have 
an immediate improvement of the maximal value of the Hilbert-Samuel 
function of the strict transform, whereas the maximal order of the weaJj 
transform stays 2. 

6.1 V{z^ + x^y^,w--> + + v'^y'^) d 

In this case the ideal to be considered is 

J/2| 33 5. 5. 32\ 
= (z + X y ,w + X + V y ) 

which is already a standard basis w.r.t. a local degree reverse lexicograph- 
ical ordering on Mon{x, y, z, w, v). 

The first auxilliary ideal A(/) — {z, x'^y^ ,x^y'^ ,w'^ , x''' ,v'^y,v'^y^) is obvi- 
ously of order 1 and dimc{^{I)/{x,y,z,w,v)^) = 1. We thus choose the 
first hypersurface of maximal contact to be defined by j/i := z. As this is 
already one of our coordinates, we do not need any coordinate change at 
this point. 

Following through the algorithmic steps of section 4, the next interesting 
degree is 5, where our corresponding ideal J has the structure: 

7 / 5 I 3 3 3 5 I 5 I 3 2i 

J — (z + z x y ,w + X + V y ), 

which is of order 5 and allows the hypersurfaces of maximal contact V{x), 
V{y), V{z), V{w), V{y) already implying that the upcoming center should 
be the origin. 

We now consider the situation after the blowing up at the origin in the 
various charts which we label (for convenience of the reader) by the gen- 
erator of the exceptional divisor on this chart: 

Chart 1: E ^ V{x) 

Istrict = (z^ + x'^y'^ + 1 + v^y"^) which can easily be checked to 
allow at most order 1 at all points of this chart. In particular, the 
order and hence the Hilbert-Samuel function have decreased under 
this blowing up. The upcoming center will be determined inside 
the hypersurface V{w^ + 1 + v'^'y^). (The weak transform of J is 
{z^ + z^x^y^,w^ + 1 + v^y^) which shows that its order has also 
dropped as was to be expected.) 

For comparison, we now observe that I-weak = {z"^ + x^y^,x^w^ + 
x^ + x'^v'^y^) where the maximal order is still 2 and the order of the 
second generator has not decreased as much as in the case of the 
strict transform. Leaving us with the first hypersurface of maximal 
contact being V{z) and a coefficient ideal of order 3 arising in this 
step. 

Chart 2: E = V{y) 

Istrict = {z'^ + x^y'^jW^ + x^ + v^) which is still of order 2 at the 
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zero locus of the ideal A(/) = {z, x^y*, x^y^ ,w'^ , ,v'^) that is along 
the line V{x, z, w, v). But the drop of order of the second generator 
causes a decrease of the maximal value of the Hilbert-Samuel, which 
is now (1, 5, 14, 29, ... ) along the line V{x, z, v, w) as compared to 
(1, 5, 14, 30, . . . ) at the previous origin. The weaJi transform of J 
is {z'' + z'^x'^y"'' ,w^' + X'' + v^) which shows that its order has also 
dropped as was to be expected. 
For comparison, we observe here that Iweak 

y^x^ + v^y"^) where the maximal order is still 2, but the order of the 
second generator has even increased producing a coefficient ideal of 
order 7 (before splitting into monomial part y^ and a non-monomial 

part of order 4) . 

Chart 3: £ = V{z) 

Istrict = (1 + x^y^z*', vet' + x^ + v^y^) of order at most 1 (analogous 
to the first chart), /^eafc = (1 + x^y^z"^, w^z^ + x^z^ + v^y^z'^) also 
of order 1. 

Chart 4: E = Viw) 

Istrict = + x^y^w'*, 1 + x^ + v^y^) of order at most 1 (again 
analogous to the first chart). Iweak = {z^ + x^y^w'^,w^ + x^w^ + 
v^y^w^) of order 2 giving rise to a coefficient ideal of order 3. 

Chart 5 E = V{v) 

Istrict = {z^ + i^t/^w'*, + a;^ + y^) of order 2, with a drop in the 
maximal value of the Hilbert-Samuel function to (1, 5, 13, 25, . . . ) 
(analogous to the second chart). Iwcak = {z^ +x'^y^v* ,w^v^ +x^v^ + 
y^v^) of order 2 giving rise to a coefficient ideal of order 5 (before 
splitting into a monomial part and a non-monomial part of order 
2). 

In this example, the order of the ideal was only two, to keep all computa- 
tions at a level of complexity which can still be followed without difficulty. 
This also made sure that forming the usual coefficient ideals of the weak 
transforms the highest powers of ideals which needed to be computed were 
second powers. 

If, however, the order of the original ideal is higher, the appearance of a 
factorial of the previous order in the exponents easily leads to far higher 
powers in the construction of the coefficient ideal which is of course iter- 
ated several times in Villamayor's approach. In our proposed approach we 
make sure that we descend the maximal possible number of hypersurfaces 
of maximal order at the very first descend of ambient dimension, thus 
minimizing the effect of taking iterated factorials. 

6.2 V{x^ + y^^,z^ + x^) (Z 

As monomial ordering we choose a negative degree reverse lexicographical 
ordering with z > y > x implying that the given generators already form a 
standard basis. Now the obvious choice of center is V{x, y, z). Considering 
the strict and weak transforms as before, we now focus on the charts in 
which the exceptional divisors are E = V{y) and E = V{z) respectively, 
omitting the discussion of the third chart: 
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Chart 1: E = V{z) 

htrict = {x^ + y^^z^', 1 + x^) is already non-singular, I^eak = {x'' + 
z^-l-a;^^*) leads to a coefficient ideal of order 8 (before splitting 
into monomial and non-monomial part) requiring up to 8-tli powers 
of certain ideals in the computation. 

Chart 2: E = V{y) 

I strict = {x^ + y^,-?^ + x^) leads by our algorithm to a locus of 
maximal Hilbert-Samuel function V{x, y, z) and will be resolved after 
the subsequent blowing up; Iweak ~ {x" + y*', y'^z'^ + y'^x^) gives rise 
to a coefficient ideal of order 10 requiring up to 60-th powers of 
certain ideals in the computation. 

Remark 7 As can he seen from the two previous examples, improvements 
of the choice of centers arise from the new variant of the algorithm, when- 
ever the lowest order generator of the given ideal is harder to resolve than 
another one, which is of significantly higher order, but 'simpler' struc- 
ture. Thus a heuristic choosing between the two approaches could take 
into account: 

• number of generators of the ideal 

• degrees of generators of the ideal 

• number of terms in lowest order generators. 

In the presence of two CPUs, another way of joining the two approaches 
could be to start them parallely each on one CPU and interrupt the slower 
one as soon as the faster one returned a result. These lines of thought, 
however, are beyond the scope of this short article and have not been ex- 
plored systematically up to now. 
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